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Navigation  Software  for  the  MPL  Vertical  Line  Array 

B.  J.  Sotirin  and  W.  S.  Hod  glass 

Marine  Physical  Laboratory 
Scripps  Institution  of  Oceanography 
La  Jolla,  CA  92093 

Introduction 

A  high  frequency  acoustic  navigation  system  is  an  integral  part  of  the  large  aperture  vertical  array 
deployed  in  4700  m  of  water  in  the  NE  Pacific  during  September  1987.  Navigation  is  defined  in  this  con¬ 
text  as  the  process  of  locating  individual  elements  of  the  array  in  3-dimension  space  at  any  particular  time. 
The  method  implemented  in  the  array  involves  a  transceiver  near  the  ocean  surface  which  sends  unique 
interrogation  signals  that  are  detected  by  bottom  moored  transponders  who  reply  with  a  pulse  at  12  KHz. 
This  reply  is  monitored  by  the  array  and  the  time  delay  between  the  initiation  and  reception  of  each  of  the 
pulses  is  calculated.  Knowing  the  sound  speed  in  the  water  column,  and  the  location  of  the  transceiver  and 
the  transponders,  this  travel  time  defines  a  slant  range  between  each  array  receiver  and  each  bottom  tran¬ 
sponder.  A  slant  range  from  a  known  point  (a  transponder)  describes  an  sphere  of  possible  receiver  loca¬ 
tions.  Intersecting  spheres  described  by  the  slant  ranges  from  three  known  points  identify  a  single  location 
if  there  are  no  errors.  Because  there  are  always  sources  of  error  in  real  data,  a  least  squares  filter  is  imple¬ 
mented  to  approximate  the  location  by  minimizing  the  squared  difference  between  the  calculated  and 
measured  values,  which  defines  the  error  in  the  assumed  position.  This  report  documents  the  data  collec¬ 
tion  methods  and  the  navigation  software  used  to  locate  12  array  receivers  during  the  September  experi¬ 
ment.  The  array  itself  is  documented  in  [Sotirin  et  al,  1988]  and  [Sotirin  and  Hildebrand,  1988].  The  array 
navigation  system,  least  squares  filter  and  navigation  data  analysis  are  documented  in  [Sotirin  and  Hilde¬ 
brand,  1989]. 

The  navigation  processing  is  separated  into  three  programs.  The  first  program  locates  the  bottom 
transponders.  The  second  program  calculates  travel  times  from  a  continuous  recording  of  the  12  KHz 
detected  returns  for  each  navigation  receiver.  The  third  program  uses  the  travel  times  and  transponder  loca¬ 
tions  output  from  the  first  two  programs  to  estimate  the  spatial  positions  of  the  receivers  using  a  non-linear 
least  squares  filter.  The  logic  flow  and  input/output  files  are  detailed  for  each  program  in  the  Appendix 
referenced  in  the  text  followed  by  a  sample  run.  The  software,  written  in  either  C  or  Fortran,  is  listed  in 
Appendix  E. 


I.  Transponder  Localization. 

The  transponder  positions  must  be  surveyed  to  acquire  the  location  parameters  defining  the  tran¬ 
sponder  net  which  are  ultimately  used  to  navigate  the  array.  Due  to  errors  in  the  measured  data,  an  estima¬ 
tion  technique  (least  squares)  must  be  implemented.  Due  to  the  nonlinear  conversion  from  travel  time 
space  to  xyz  positional  space,  the  least  squares  method  proceeds  iteratively.  A  data  set  is  obtained  contain¬ 
ing  spatial  positions  and  travel  time  measurements.  If  all  parameters  were  known  exactly,  the  travel  times 


calculated  using  die  spatial  positions  and  the  measured  value  would  be  identical.  This  is  obviously  not  the 
case,  and  the  difference  between  the  calculated  and  measured  values  defines  the  error  in  the  assumed  posi¬ 
tion  which  is  minimized  during  the  iteration.  Transponder  positions  with  accuracies  of  less  than  a  meter 
are  achieved  by  this  method.  The  description  below  details  the  data  collection,  the  program  inputs,  the 
least  squares  implementation,  and  the  resulting  transponder  positions. 

Data  Collection.  The  input  data  set  for  the  least  squares  iteration  requires  initial  spatial  positions  of 
the  transponders  and  transceiver  in  meters  from  an  arbitrary  origin,  and  slant  ranges  between  the  tran¬ 
sponders  and  transceiver.  The  data  collected  in  the  form  of  travel  times  and  Global  Positioning  Satellite 
(GPS)  fixes  are  transformed  into  the  input  parameters  required  by  the  least  squares  filter.  GPS  fixes  are 
converted  into  initial  xy  positions  in  meters  from  an  arbitrary  origin.  Travel  times  are  converted  to  slant 
ranges  with  knowledge  of  the  local  sound  speed  profile.  The  travel  times  required  are  normally  recorded 
during  an  intensive  surface  ship  survey  during  which  the  ship  criss-crosses  the  area  in  which  the  tran¬ 
sponders  were  deployed  recording  the  travel  time  data  and  GPS  fixes  for  its  own  position.  [Spiess,  1985], 
[Smith  et  al,  1975]  A  transceiver  is  either  hull  mounted  or  towed  on  a  short  line  such  that  its  position  rela¬ 
tive  to  the  ship  is  known.  The  transceiver  sends  a  continuous  stream  of  unique  transponder  interrogation 
pulses,  and  the  transponder  replies  are  recorded  while  the  ship  criss-crosses  the  area.  Using  this  method,  a 
series  of  travel  times  are  obtained  from  a  wide  variety  of  ship  positions.  As  the  ship  crosses  over  the  top  of 
a  transponder,  an  estimate  of  transponder  depth  is  acquired,  and  as  the  transponder  baselines  are  crossed, 
intertransponder  distances  are  defined. 

During  the  September  experiment,  although  travel  times  were  recorded  as  described,  the  12  KHz 
receiver  was  deployed  on  a  200  m  line  due  to  the  noise  level  of  the  ship.  This  introduces  the 
interrogator/receiver  depth  as  an  unknown  defining  an  underdetermined  set  of  equations  which  cannot  be 
solved  for  a  unique  solution.  Since  the  measurements  described  above  were  not  sufficient  for  the  tran¬ 
sponder  survey,  the  travel  time  data  collected  by  the  navigation  equipment  installed  on  FLIP  were  utilized 
instead.  This  is  unfortunately  not  an  optimum  choice  of  observation  configurations  [Spiess,  1985]  so  the 
results  of  several  simulations  devised  to  assess  the  effect  on  the  estimated  positions  are  also  presented.  The 
horizontal  motion  of  FLIP  was  constrained  by  a  three  point  moor.  The  FLIP  data  set  provides  sufficient 
range  information  but  the  azimuthal  component  was  not  well  constrained.  The  vertical  component  was 
estimated  from  the  echo  sounding  depth  at  FLIP  and  apriori  information  that  the  sea  floor  in  the  experiment 
area  was  relatively  flat  The  inputs  used  to  locate  the  transponders  are  the  vertical  sound  speed  profile  in 
the  test  area,  the  initial  xy  positions  of  FLIP  from  the  arbitrary  origin,  the  depth  of  the  interrogation  tran¬ 
sponder  hardwired  o  the  FLIP ,  the  slant  ranges  from  each  transponder  to  FLIP  as  determined  by  chart 
recorder  traces  and  associated  errors,  the  initial  transponder  xy  positions  from  the  arbitrary  origin  using  the 
satellite  fixes  and  the  transponder  depths. 

The  initial  GPS  positions  in  latitude  and  longitude  are  converted  to  xy  distances  in  meters  from  an 
arbitrary  origin.  For  the  September  experiment,  the  origin  was  chosen  as  34°  47’  N,  126°  00’  W,  with  x 
increasing  positively  toward  the  east  and  y  increasing  positively  toward  the  north.  A  GPS  position  of  the 
ship  is  recorded  during  each  transponder  deployment  The  transponders  have  45  kg  negative  buoyancy 
when  they  are  deployed  from  the  fantail  of  the  surface  ship  to  insure  that  the  GPS  position  of  the  ship  is  an 
accurate  initial  position  of  the  actual  transponder  position  on  the  sea  floor.  Latitude  and  longitude  for  each 
set  of  travel  time  measurements  (one  measurement  from  each  transponder)  are  also  converted  to  m  to  pro¬ 
vide  the  initial  estimate  of  the  transceiver  position  during  the  survey.  The  ’survey’  from  FLIP  had  such  a 
limited  range  that  the  same  initial  position  was  used  for  all  measurement  sets.  During  a  normal  ship  sur¬ 
vey,  travel  time  data  is  collected  from  a  large  variety  of  horizontal  positions  and  individual  initial  positions 
are  important.  If  GPS  positions  are  not  available,  a  technique  detailed  in  Appendix  D  may  be  used  to  esti¬ 
mate  the  initial  position  from  the  measured  travel  time  data.  Positions  in  latitude  and  longitude  were  con- 
voted  to  meters  by  calculating  the  radius  of  the  earth  at  the  position  using  the  following  equations  [Sta¬ 
cey]: 
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x  =  rp  &{lat ) ,  y  =rp  cos  (lat )  6 (long ) 
rp=a(l-f  sinKlat)),  f= 

where  x  is  the  E-W  distance  in  meters  of  the  position  from  the  origin,  y  is  the  N-S  distance  in  meters  of  the 
position  from  the  origin,  rp  is  the  radius  of  the  earth  at  the  position  latitude,  Ulat )  is  the  difference  in  lati¬ 
tude  in  radians  between  the  position  and  the  origin,  8 {long )  is  similarly  the  difference  in  longitude,  lat  is 
the  latitude  of  the  position,  a  is  the  equitorial  radius  of  the  earth  in  meters  and  c  is  similarly  the  polar 
radius. 

The  conversion  between  travel  time  and  xyz  positions  requires  knowledge  of  the  sound  speed  profile 
at  the  experiment  site.  The  local  sound  speed  profile  was  calculated  from  measurements  of  conductivity, 
temperature  and  depth.  Travel  time  deviations  due  to  the  variations  within  the  thermocline  and  to  refrac¬ 
tion  of  acoustical  energy  were  shown  to  be  negligible  for  this  experiment  therefore  a  constant  sound  speed 
(harmonic  mean)  was  used  for  the  conversion.  [Sotirin  and  Hildebrand,  1989] 

The  travel  time  measurements  used  for  transponder  localization  were  collected  by  transmitting  the 
transponder  interrogation  signals  from  FLIP  once  an  hour  for  18  days  and  recording  the  returns  on  a  chart 
recorder.  The  transceiver  is  mounted  on  the  bottom  (90  m  in  depth)  of  FLIP  and  the  pulse  level  is  adjusted 
manually  above  the  ambient  noise  for  consistent  transponder  replies.  Round  trip  travel  times  for  navigating 
FLIP  are  measured  carefully  by  hand  on  the  chart  recorder  output  with  an  estimated  random  gaussian  error 
of  2-3  ms.  The  chart  recorder  trace,  set  on  a  one  second  sweep  rate,  records  the  filtered  12  kHz  (500  Hz 
bandwidth)  replies  received  from  the  transponder  being  interrogated  (Figure  1.1)  delayed  6-7  seconds  from 
the  interrogation  pulse.  Each  transponder  is  interrogated  individually  by  transmitting  its  unique  signal,  trig¬ 
gered  by  the  chart  recorder,  once  per  second  for  45  seconds,  and  notating  its  reply  on  the  chart  by  sweep 
number  with  color-coded  pens  (red,  green  or  blue).  The  figure  shows  the  direct  and  multipath 
(surface/bottom  bounce)  returns  for  each  transponder.  Differentiation  between  the  surface  and  bottom 
bounces  is  difficult  because  the  FLIP  transceiver  is  the  same  distance  below  the  surface  as  the  transponder 
is  above  the  sea  floor.  The  direct  return  from  the  red  and  blue  transponders  is  strong  and  consistent. 

Several  error  modes  are  evident  in  the  returns  from  the  green  transponder  however.  The  sporatic  direct 
return  of  the  green  transponder  (Gj)  is  caused  by  the  transponder  detecting  either  the  surface  or  bottom 
bounce  of  the  transmitted  signal.  The  error  was  attributed  to  the  transponder  because  the  same  transceiver 
system  detects  the  returns  from  all  transponders  and  this  type  of  error  was  not  evident  on  all  transponders. 
Detection  of  the  multipath  arrival  of  the  interrogation  pulse  delays  the  arrival  of  the  transponder  direct 
return  such  that  it  arrives  at  the  same  time  as  the  multipath  return  of  a  direct  interrogation  pulse  (G„i).  The 
multipath  return  of  the  multipath  interrogation  arrives  120  ms  later  {Gmi)  clearly  identifying  the  first  error 
mode.  The  second  error  mode  of  the  green  transponder  is  displayed  as  the  smaller  perturbations  in  the  time 
of  arrival  of  the  Gd  return.  It  is  more  difficult  to  identify  because  the  deviation  from  the  true  return  is 
smaller  (although  still  significant)  and  not  as  consistent  as  the  first  type  of  error.  These  error  modes  were 
identified  during  about  20%  of  the  experiment  at  random  intervals,  causing  concern  during  array  detection 
post-processing  discussed  in  the  next  section. 

Software  Implementation.  Once  the  ingredients  for  the  least  squares  method  have  been  accumu¬ 
lated,  all  xyz  positions  are  adjusted  until  the  root  mean  squared  (rms)  error  satisfies  the  convergence  cri¬ 
teria.  This  is  accomplished  in  several  parts  by  first  maintaining  constant  transponder  positions  and  perturb¬ 
ing  the  FLIP  positions,  then  holding  the  current  FLIP  positions  constant  while  perturbing  the  transponder 
positions,  and  finally  examining  the  mean  squared  error  of  each  FLIP  position  to  determine  whether  it 
should  be  preserved  as  a  viable  contributor.  A  general  outline  of  the  method  is  presented  in  Figure  1.2  and 
discussed  below.  This  is  followed  by  a  detailed  description  of  the  conversion  from  travel  times  to  slant 
ranges  and  spatial  positions. 

The  program  consists  of  three  concentric  stages,  the  first  stage  adjusts  the  xy  positions  in  two  inner 
loops  each  testing  the  rms  error  against  specific  convergence  criteria  for  the  least  squares  filter,  the  second 
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stage  tests  the  total  transponder  rms  error  against  a  second  set  of  convergence  criteria,  and  the  third  stage 
tests  yet  again.  The  three  stagesare  referred  to  within  the  software  as  stage  1:  xpfil,  stage  2:  xploop  and 
stage  3:  xpmcdn.  The  two  first  stage  loops  shown  in  Figure  1.2  have  similar  internal  operations.  The  meas¬ 
ured  travel  times,  sound  speed  profile,  FLIP  depth  and  transponder  depths  are  considered  known,  while  the 
xy  positions  of  FLIP  and  the  transponders  are  considered  unknown.  The  measured  travel  times  are  con¬ 
verted  to  slant  ranges  by  multiplying  by  the  harmonic  mean  of  the  sound  speed  profile.  This  slant  range  is 
projected  as  a  horizontal  range  and  compared  to  the  range  calculated  from  the  xy  positions.  For  the  first 
loop,  the  squared  difference  between  these  two  ranges  are  summed  over  the  number  of  transponders  for  a 
particular  FLIP  position  and  the  square  root  is  taken  to  yield  the  rms  error.  If  the  rms  error  does  not  satisfy 
the  convergence  criteria,  then  the  FLIP  position  is  adjusted  and  the  loop  repeats  with  the  adjusted  position. 
The  adjustment  is  calculated  with  the  search  direction  as  the  negadve  gradient  and  the  step  size  as  a  con¬ 
stant  (1.5  m)  unless  tire  rms  error  is  less  than  1  m  at  which  time  the  step  size  begins  to  decrease;  as  the 
minimum  is  approached,  the  step  size  is  calculated  as  a  function  of  the  percent  change  in  interned  rms 
error.  The  convergence  criteria  for  the  first  stage  loops  are  defined  such  that  an  absolute  rms  error  less  than 
0.15  m,  a  0.015  percent  change  in  the  iterated  rms  error,  or  a  maximum  number  of  iterations  (30)  will  ter¬ 
minate  the  loop  and  save  the  current  position.  These  criteria  test  each  adjusted  position  individually.  The 
first  loop  is  repeated  for  each  FLIP  position.  The  second  loop  of  the  first  stage  performs  the  same  cadence 
maintaining  the  current  FLIP  positions  constant  while  adjusting  the  transponder  positions.  The  rms  error  is 
calculated  over  the  number  of  FLIP  positions,  and  evaluated  using  the  same  criteria.  The  second  loop  is 
repeated  for  each  transponder  position.  In  the  second  stage,  these  two  loops  are  initiated  again  based  on 
the  percent  reduction  (<  0.35%)  in  the  transponder  rms  error  summed  over  all  transponders  and  total 
number  of  iterations  (>  30).  Upon  completion  of  the  first  two  stages,  the  rms  error  for  all  the  transponders 
is  calculated  and  evaluated  according  to  the  following  criteria: 

1)  The  rms  error  *  ^number  of  FLIP  positions  is  <  1.0. 

2)  The  percent  reduction  in  rms  error  is  less  than  a  specified  value  (0. 1  %). 

3)  The  absolute  rms  error  is  less  than  a  specified  number  (0.75). 

If  the  rms  error  satisfies  any  of  the  above  criteria,  the  current  positions  are  written  out  and  the  program  is 
terminated.  If  the  rms  error  does  not  satisfy  any  of  the  above  criteria,  then  the  rms  errors  associated  with 
each  FLIP  position  are  examined  and  any  position  with  an  error  greater  than  the  transponder  rms  error 
times  a  user  specified  value  (default=2)  is  deleted  from  the  array  and  the  entire  process  begins  again.  The 
rms  error  is  initialized  to  10000  prior  to  each  loop/stage  so  that  unless  the  absolute  error  is  small,  the 
loop/stage  will  always  be  executed  more  than  once. 

The  conversion  from  travel  time  to  slant  range  and  the  spatial  position  adjustment  is  straight  forward 
and  executed  within  the  subroutines  xxcor  and  xpfil.  A  detailed  outline  of  the  program  flow  is  found  in 
Appendix  A.  The  program  originates  with  the  Marine  Physical  Laboratory’s  DEEPTOW  group  and  has 
been  in  existance  since  the  late  1960’s.  The  version  documented  here  is  the  most  recent  in  a  long  series 
and,  unfortunately,  every  programmer  has  left  a  marie.  Consequently,  the  subroutines  each  have  unique 
variable  names  for  the  parameters,  increasing  the  confusion  in  documentation.  The  text  below  attempts  to 
maintain  the  variable  names  within  the  subroutines  specified.  The  travel  time  to  slant  range  conversion  is 
computed  in  xxcor.  The  measured  information  is  input  to  the  program  as  a  slant  range  assuming  a  homo¬ 
geneous  medium  with  a  sound  speed  of  1500  m/s.  This  is  not  normally  the  case,  and  for  the  September  test 
the  sound  speed  profile  was  measured  and  input  as  a  data  statement  in  xxcor  into  the  array  vdp  as  a  hor¬ 
izontally  stratified  medium.  If  sound  speed  corrections  are  requested  by  the  user,  the  input  slant  ranges  are 
converted  back  to  the  original  time  measurement  (t)  and  the  slant  range  contribution  for  each  sound  speed 
layer  is  summed,  returning  the  corrected  slant  range  (xnew).  If  the  slant  range  is  a  depth,  for  example,  the 
summation  starts  at  the  surface  and  sums  the  contribution  in  each  layer  until  the  accumulated  time  (r„)  is 
equal  to  the  measured  travel  time.  This  is  an  implementation  of  the  following  equation  in  which  x  is  the 
corrected  slant  range: 
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Figure  1.2  Transponder  Localization  Program  Flow.  A  general  outline  of  3  stages  in  the  least  squares  program 
logic  is  presented.  During  stage  1,  the  FLIP  positions  are  iterated  first,  followed  by  the  transponder  positions  and 
the  rms  errors  are  calculated  for  the  positions  iterated  and  compared  to  a  set  of  convergence  criteria;  stage  2  com¬ 
pares  transponder  rms  errors  to  a  second  set  of  convergence  criteria;  stage  3  compares  transponder  rms  errors 
against  a  third  set  of  convergence  criteria  and  determines  whether  to  delete  FLIP  positions  with  excessive  errors. 
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where  in  terms  of  the  subroutine,  20  =  xdepth  =  0, 27  =  z  is  the  unknown  depth,  C(z)  =  v  is  the  sound  speed 
in  the  layer,  and  dz  =  deltad  is  the  layer  depth.  When  the  unknown  is  a  slant  range,  the  implementation  is 
an  approximation  to  the  above  equation  (where  dz  is  the  slant  range  component  in  the  layer  assuming  the 
launch  angle  is  constant)  which  is  accurate  for  the  configuration  of  the  data  set  considered  for  the  Sep¬ 
tember  experiment  Errors  increase  as  the  horizontal  range  or  projection  (defined  below)  of  the  slant  range 
increases  however,  so  for  larger  scale  experiments  a  more  accurate  implementation  is  advised. 

Once  the  slant  range  has  been  calculated,  the  conversion  to  spatial  coordinates  is  visualized  as  the 
geometric  relationship  between  the  transmitter  (t  subscript)  and  receiver  (r  subscript): 

slant  range  =  [  ( x ,  -  )2  +  (y,  -  yr)2  +  (2,  -  zr )2  ]T. 

Because  the  error  in  the  depth  parameters  are  small  compared  to  the  horizontal  parameters,  the  slant  range 
is  projected  onto  the  xy  plane  prior  to  the  adjustment  (Figure  1.3): 

1  1 

horizontal  projection  =  [(slant  range  )2  -  (z,  -  z,)2]7  =  [(z,  -xr)1  +  (y,  -  yr)2]7. 

1 

The  left  half  of  the  equation  is  calculated  in  xpread,  CRANS (NTR  NPOS)  =  (S2-/)2)7  where  NTR  =  the 
number  of  transponders,  NPOS  -  the  number  of  FLIP  positions,  S  =  slant  range  between  a  position  and  a 
transponder,  D  =  transponder  depth  -  FLIP  depth,  and  passed  into  xpfil  as  an  array  IIRAN  which  is 
redefined  within  a  loop  as  a  variable  HH .  The  right  half  of  the  equation  is  calculated  in  xpfil, 

RNGEC  =  ^(XDIFF )2  +  (YDIFF )2,  where  XDIFF  =XG  -  XF  (NDAT),  YDIFF  -YG  -  YF (NDAT), 

(XG ,  YG )  are  the  xy  positions  being  interated,  (XF ,  YF )  are  the  fixed  positions  indexed  over  NDAT  A  data 
points,  and  the  error  ERR ,  summed  over  the  number  of  fixed  positions,  is  (RNGEC -HH)2. 

The  adjustment  is  calculated  using  the  steepest  descent  method  to  minimize  the  error  by  following 
the  mean  squared  error  gradient  to  a  minimum.  For  known  transponder  positions  and  the  x-direction,  the 
perturbed  position  is: 

ndata  _ 

XG  =  _v  XG  +  h  ERR  (NDAT) 

NuXi  “1 

where  h  is  the  step  size,  and  ERR’  is  the  negative  derivative  of  the  error  function  with  respect  to  XG.  The 
y-direction  adjustment  is  calculated  similarly.  The  error  derivative  expands  to: 

ERR’=  d(RNGEC-HH)2  m 2(////  _ MGEc)MPlFF2+YDlFFy 


_  (HH  -RNGEC)  d(XG-XF)2 
RNGEC —  JXG 


=RATIO*XDIFF 


where  RATIO  =  (HH  -  RNGEC  )/RNGEC ,  and  the  constants  are  absorbed  by  the  step  size  h. 


Simulations.  Several  simulations  of  the  array  navigation  system  were  conducted  to  examine  the  sen¬ 
sitivity  of  the  estimated  transponder  positions  to  errors  in  travel  time  measurements  and  initial  positions. 
The  spatial  configuration  used  closely  resembles  the  experimental  set  up  of  the  September  sea  test  as 
shown  in  Figure  1.4.  The  transponders  and  FLIP  were  initially  assigned  to  known  positions  with  deter¬ 
mined  slant  ranges  as  shown.  Two  simulations  were  conducted  to  illustrate  the  transponder  position 
response  to  errors;  two  other  simulations  were  conducted  to  show  the  effect  of  transponder  and  FLIP  posi¬ 
tional  errors  on  the  estimated  array  positions  and  are  presented  in  [Sotirin  and  Hildebrand,  1989].  The 
result  of  the  first  simulation  was  the  3D  error  surface  for  various  parameters.  The  second  was  a  Monte 
Carlo  simulation  of  the  initial  transponder  positions.  The  simulations  provide  an  understanding  of  the 
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Figure  1.3  Navigation  Overview.  The  horizontal  projection  is  estimated  first  by  using  the  measured  slant  range 
and  depths  ( HRAN)  and  then  by  using  the  initial  xy  positions  (mgxy).  The  initial  xy  positions  of  the  transponders 
(Tl=red,  T2=green,  T3=blue)  and  FLIP  as  measured  by  a  GPS  fix  are  notated  in  meters  on  the  axis,  and  beneath  the 
transponder  for  depth. 
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Figure  1.4  Spatial  Configuration  for  Navigation  During  the  September  1987  Sea  Test.  The  estimated  xy  posi¬ 
tions  of  the  3  fixed  bottom  transponders  and  of  FLIP  measured  hourly  over  18  days  are  shown  in  plan  view.  The 
mooring  lines  are  represented  by  arrows,  and  the  slant  range  (s)  and  horizontal  projection  of  the  slant  range  (h)  are 
indicated  in  meters  from  an  arbitrary  FUP  xy  position  of  (24003200)  m.  Transponder  baseline  distances  are  also 
indicated  in  meters. 
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effect  of  the  unconventional  survey  data. 

To  examine  the  error  surface  for  the  least  squares  configuration  used  in  the  previous  section,  specific 
model  parameters  were  perturbed  systematically.  The  error  is  the  rms  value  of  the  differences  between  the 
slant  ranges  based  on  the  travel  time  measurements  corrected  for  sound  speed  deviations  (Eq.  1.1  dividing 
through  by  c)  and  those  calculated  from  the  spatial  positions  of  the  transponders  and  FLIP  which  were  out¬ 
put  from  the  least  squares  filter.  Ideally  this  produces  a  single  well-defined  minimum  for  which  the  optimi¬ 
zation  method  searches.  This  is  not  the  case  in  many  real  applications  however,  and  the  possibility  of  local 
minima  and/or  a  broad  global  minimum  should  be  investigated.  The  inverted  error  surface  for  perturba¬ 
tions  in  the  horizontal  positions  of  the  blue  transponder  is  shown  in  Figure  1.5.  It  exhibits  a  narrow  chan¬ 
nel  of  local  minima  which  appear  as  a  ridge  plotted  as  the  negative  logarithm  of  the  error  against  the  per¬ 
turbation  amplitude  in  x  and  y  position.  The  ridge  is  orientated  perpendicularly  to  the  FL/F/transponder 
range  direction  indicating  that  the  azimuthal  component  is  not  well  constrained.  Perturbing  the  horizontal 
positions  of  the  other  transponders  produced  similar  results. 

The  system  displays  a  sensitivity  to  initial  positions  which  is  shown  to  be  an  artifact  of  the  limitations 
in  the  FLIP  ’survey’  discussed  previously  and  of  the  structure  of  the  error  surface.  The  direction  in  which 
the  transponders  are  moved  is  constrained  to  nearly  parallel  to  the  FL/F-transponder  baseline  (Figure  1.6). 
This  occurrs  because  the  distribution  of  FLIP  positions  shown  in  Figure  1.4,  provide  adequate  range  infor¬ 
mation  but  minimal  azimuthal  information  regarding  the  transponder  positions  as  was  indicated  by  the 
error  surface.  There  are  no  sure  techniques  for  locating  a  global  minimum  in  the  company  of  local  minim- 
ina.  Without  independent  positional  data  corroborating  the  results,  either  Monte  Carlo  techniques  would 
have  to  be  incorporated  in  the  initial  position  of  each  transponder,  or  knowledge  of  the  error  surface  would 
have  to  be  employed  as  a  constraint  (search  range/azimuth  space  rather  than  xy  space).  Fortunately,  the 
initial  estimate  of  the  transponder  positions  were  acquired  from  accurate  GPS  fixes  and  the  resulting  error 
in  GPS  positions  compared  to  the  estimated  FUP  positions  at  corresponding  times  had  an  rms  value  of 
only  10  m.  Thus  the  estimated  transponder  positions  were  declared  adequate. 


II.  Array  Travel  Time  Acquisition. 

Travel  time  measurements  were  acquired  by  interrogating  three  bottom  mounted  transponders  from 
FLIP  and  detecting  their  replies  at  the  navigation  receivers  distributed  across  the  900  m  aperture  array. 
There  are  24  navigation  receivers  located  at  ±  3.75  m  from  each  processor  which  are  separated  by  75  m 
(Figure  2.1).  Due  to  bandwidth  constraints,  data  from  12  of  the  navigation  receivers  (one/section)  were 
decimated  and  recorded  during  the  September  1987  sea  test.  The  data  bit  stream  from  transmit  time  of  the 
interrogation  pulse  to  receive  time  of  the  reply  at  the  array  was  reconstructed  during  post-processing  for 
each  navigation  receiver  and  the  travel  times  calculated. 

The  navigation  timing  was  based  on  a  16  bit  clock  driven  at  a  1  KHz  rate  which  initiated  the  naviga¬ 
tion  sequence.  This  clock  is  referred  to  as  the  hardware  clock  to  differentiate  it  from  the  real-time  clock 
(local  time)  and  the  Greenwich  Mean  Time  (GMT)  clock.  The  timing  is  illustrated  in  Figure  2.2  and 
described  below.  A  transceiver  located  at  the  bottom  of  FLIP  transmited  a  series  of  65536  ms  transponder 
sequences  (Figure  2.2a).  A  transponder  sequence  consisted  of  4  interrogation  pulses  10  ms  long  at  10  s 
intervals  beginning  at  a  16  bit  clock  rollover  followed  by  a  35.536  s  silent  interval  (Figure  2.2b).  The  first 
three  pulses  were  at  the  unique  transponder  interrogate  frequencies  of  the  bottom  transponders  (10, 10.5 
and  1 1  KHz).  Upon  receiving  an  interrogation  pulse  the  bottom  transponders  replied  with  a  3  ms  pulse  at 
12  KHz.  The  turn  around  time  of  the  transponders  is  signal  to  noise  dependent,  however  due  to  the 
strength  of  the  transmitted  signal,  the  delay  was  assumed  to  be  less  than  1  ms.  The  fourth  interrogation 
pulse  transmitted  by  the  FLIP  interrogate  transceiver  was  at  12  KHz  which  simulates  a  bottom  transponder 
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Figure  1 .5  3D  Error  Surface.  The  estimated  positions  of  FLIP  and  the  transponders  were  considered  known, 
selected  parameters  were  perturbed  and  the  resulting  error  calculated.  In  this  example,  the  blue  transponder  x  and  y 
positions  were  perturbed  in  1  meter  increments  to  ±100  m  from  the  original  known  position  which  is  plotted  at  the 
center.  The  z  axis  is  plotted  as  the  negative  logarithm  of  the  error  to  allow  visualization  of  the  minima. 


11 


METERS  EAST  FROM  126°  O'  W 


Figure  1.6  Monte  Carlo  simulation  of  initial  transponder  positions.  This  simulation  illustrates  the  constraints  in 
search  direction  placed  on  the  least  squares  iteration  by  the  spatial  configuration  of  the  survey  points.  The  azimuthal 
position  of  the  transponders  is  not  accurately  determined  by  the  FLIP  data  set  used.  The  transponders  tend  to  move 
along  the  transponder  to  FLIP  direction.  The  initial  transonder  position  is  notated  by  a  triangle,  the  final  position  by 
across. 
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Figure  2.2  Navigation  Timing  Diagram.  The  timing  associated  with  the  array  navigation  system  is  illustrated, 
(a)  The  transponder  interrogation  sequence  is  initiated  every  65536  (64  K)  ms  by  a  16-bit  hardware  clock  rollover 
(R).  (b)  Each  transponder  sequence  consists  of  4  transmissions  at  unique  frequencies  issued  from  a  transceiver 
mounted  on  the  bottom  of  FLIP  (90  m  depth)  spaced  10  s  apart  (c)  The  array  receives  the  12  kHz  transponder 
replies  and  the  12  kHz  pulse  from  FLIP.  The  time  of  arrival  of  a  transponder  reply  reflects  the  time  it  took  the 
acoustic  pulse  to  travel  from  FLIP  to  the  transponder  to  the  array  receiver.  The  time  of  arrival  of  the  FLIP  12  kHz 
pulse  corresponds  to  the  distance  from  the  array  receiver  to  FLIP,  (d)  The  hardware  clock  rollover  and  12  kHz 
arrivals  may  occur  anywhere  within  a  128  ms  tape  buffer  which  complicates  the  data  extraction  procedure. 


reply.  The  array  therefore  received  four  consecutive  12  KHz  reply  pulses  about  once  a  minute  (Figure 
2.2c).  The  navigation  receivers  were  capable  of  detecting  a  12  KHz  signal  6  dB  below  the  noise  level  in  a 
200  Hz  band.  The  binary  output  of  the  detector  was  sampled  at  5  KHz  and  decimated  by  2  to  provide  a 
continuous  time  series  consistent  with  the  5  bits  per  processor  every  2  ms  allowed  for  navigation  data 
within  the  specified  data  format  Thus  the  data  from  12  navigation  receivers  located  3.75  m  below  each 
processor  (H5  in  Figure  2.1)  in  the  array  were  multiplexed  in  with  the  low  frequency  acoustic  data, 
transmitted  to  the  surface  and  recoded.  The  interrogation  sequence  was  synchronized  with  the  timebase  in 
the  array  and  the  initiation  time  of  the  sequence  (as  indicated  by  the  rollover  of  the  hardware  clock)  was 
sampled  every  128  ms  and  recorded  on  the  tape.  The  tape  format  showing  the  placement  of  the  navigation 
data  is  illustrated  in  Figure  2.3.  The  navigation  data  time  series  are  reconstructed  from  this  recorded  data 
by  the  method  described  below. 

The  structure  of  the  data  on  tape  was  not  optimized  to  facilitate  extraction  of  the  navigation  bit 
stream.  The  header  containing  the  timing  information  appears  in  the  first  8  words  of  each  tape  buffer 
which  contains  128  ms  of  data.  Most  of  the  6280  16-bit  words  within  a  buffer  are  assigned  to  low  fre¬ 
quency  acoustics  rather  than  navigation  as  seen  in  Figure  2.3.  To  reconstruct  the  sampled  dme  series  out¬ 
put  from  each  navigation  receiver,  the  5  bits/processor  must  be  extracted  and  stored.  The  major  program¬ 
ming  effort  was  in  initializing  and  incrementing  pointers  and  in  error  checking.  A  description  of  the  main 
program  and  subroutines  and  a  sample  run  will  be  found  in  Appendix  B.  The  extraction  of  the  data  begins 
with  the  data  buffer  containing  the  hardware  clock  rollover.  The  rollover  may  occur  anywhere  within  the 
data  buffer.  The  program  examines  each  buffer  header  for  the  hardware  clock  rollover.  This  is  determined 
by  testing  between  consecutive  clock  times  for  a  negative  difference,  and  linearly  interpolating  within  the 
buffer  for  the  correct  frame.  For  example,  if  the  hardware  clock  in  buffer  1  is  hwc  i  =  65430  and  the 
hardware  clock  in  the  next  consecutive  buffer  is  kwc2  =  22,  then  the  difference  (hwc2-  hwc\  =  -65408)  is 
negative  and  the  clock  rollover  occurs  within  buffer  1  at  frame  53;  the  data  stored  would  begin  at  frame  0, 
buffer  1.  The  time  difference  from  the  beginning  of  the  buffer  to  the  rollover  is  stored  in  a  parameter 
called  start  and  passed  to  the  subroutine  navloc  which  calculates  the  travel  time.  Once  the  rollover  is 
identified,  the  navigation  data  associated  with  each  processor  is  masked  off  and  stored  as  the  10  most 
significant  bits  in  a  16-bit  word  (fillnav ). 

The  travel  time  calculation  is  a  simple  difference  once  the  time  of  the  transponder  reply  is  deter¬ 
mined.  To  locate  the  reply,  the  data  is  treated  as  a  time  series  of  0.4  ms  bits,  and  a  correlation  between  the 
time  series  data  and  a  replica  transponder  pulse  is  initiated  as  a  matched  filter  detector.  Due  to  inconsisten¬ 
cies  in  receiver  detection  threshold  and  noise  level,  each  receiver  is  assigned  an  individual  replica  pulse 
length.  Temperature  sensitivity  of  the  capacitors  and  high  failure  rate  of  the  inductors  in  a  phase  matching 
tuned  filter  caused  the  mismatch  in  the  detection  threshold  of  the  individual  receivers.  This  variation  in 
receiver  error  is  apparent  in  the  positional  error  distribution  shown  in  Figure  2.4,  in  which  the  only 
hardware  or  processing  difference  is  in  the  navigation  receiver  itself.  Additional  variation  in  receiver  sig¬ 
nal  to  noise  level  of  the  incoming  reply  could  have  been  caused  by  the  12  KHz  beam  pattern  of  the  indivi¬ 
dual  elements.  Each  element  consisted  of  two  hydrophones  wired  in  series  with  a  spacing  between  8  and  9 
cm  such  that  at  12  KHz,  a  notch  in  the  beam  pattern  appears  between  46°  and  52°  from  broadside.  From 
Figure  1.4,  the  arrival  angle  of  the  transponder  replies  when  the  array  was  nominally  at  400  m  was  between 
58°  and  65°  from  broadside  and  signal  to  noise  was  high  enough  to  be  detected;  at  deeper  depths,  however, 
this  is  potentially  a  problem.  The  rms  error  used  in  the  distribution  calculations  accompanies  each  receiver 
position  estimation  described  in  Section  III.  The  transponder  reply  is  a  CW  pulse,  so  the  correlator  is 
implemented  as  a  moving  adder,  with  a  detection  defined  as  the  first  occurrance  of  a  normalized  correlation 
amplitude  greater  than  or  equal  to  1.0  within  a  valid  data  window. 

The  window  was  installed  due  to  excessive  noise  levels  and  the  interference  seen  in  the  return  of  a 
single  navigation  receiver  from  the  hardware  clock  rollover  as  shown  in  Figure  2.5.  The  window  must  be 
determined  prior  to  program  execution.  An  option  in  the  program  harrynav  (-p)  will  print  out  all  bits  set  by 
the  navigation  receiver,  packed  into  the  least  significant  10  bits  of  a  16-bit  word  (for  a  maximum  of  3FF)  to 
allow  the  user  to  determine  the  window  parameter.  An  example  is  seen  in  Figure  2.6:  the  return  for 


15 


16 


Figure  2.3  Data  Tape  Format.  The  tape  consists  of  11000-12000  buffers,  a  buffer  includes  an  8  word  header  and  64  98-word  frames,  the  header  stores 
relevant  timing  and  error  information,  a  frame  contains  a  frame  sync  word,  a  frame  counter,  and  twelve  8-word  processor  data  groups,  a  processor  data  group  has 
ten  samples  of  the  low  frequency  acoustic  field,  5  navigation  bits  and  a  3  bit  processor  ID.  The  information  relevant  to  the  navigation  includes  the  8  word 
header,  frame  sync  word,  telemetry  module  word,  5  (0.4  ms  sample  rate)  navigation  bits  (bit  8  occuring  first  in  real  time)  and  a  3  bit  processor  ID  (PID). 


Observations /Bin 
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Figure  2.5  Correlated  Output  from  a  Single  Receiver.  A  single  receiver  should  ideally  display  only  4  detections 
during  a  transponder  interrogation  sequence.  Shown  is  a  40  s  time  series  of  the  correlated  output  of  processor  3. 
Each  line  represents  one  second  of  data,  a  signal  is  detected  as  a  transponder  return  if  the  amplitude  is  above  the 
zero  level  of  the  line  above  it  The  noise  and  interferring  signals  warren  ted  that  only  data  within  a  valid  data  win¬ 
dow  be  considered.  The  desired  detections  are  indicated  by  the  stipled  boxes. 


processor  2  is  clean  and  the  transponder  reply  is  quite  clear.  The  last  4  bits  of  word  i503  are  set  followed 
by  8  bits  in  word  606,  for  a  total  of  12  sequential  bits  set.  The  words  have  a  maximum  of  10  bits  or  4  ms. 
The  valid  data  window  length  is  a  minimum  of  100  ms  and  varies  with  the  transponder.  Only  the  leading 
edge  of  the  window  is  specified  in  the  window  parameter  input  file  window.dat.  An  estimate  of  the  leading 
edge  of  the  window  for  a  particular  receiver  p  may  be  obtained  by  multiplying  the  first  word  of  the  tran¬ 
sponder  reply  Wp  by  4  and  subtracting  184  ms.  A  window  parameter  for  each  processor  and  each  tran¬ 
sponder  must  be  specified,  48  total.  The  algorithm  for  computing  the  time  of  arrival  of  the  leading  edge  of 
the  reply  is  rfl  =Wp  *  4  +  (10  -  )  *  0.4  -  start  -c  where  Wp  is  the  first  word,  and  Bp  is  the  number  of 

bits  set  in  the  first  word,  start  is  the  time  difference  between  the  beginning  of  the  data  buffer  and  the 
hardware  clock  rollover  as  described  earlier  (printed  in  line  1  of  Figure  2.6  as  Start:  14),  and  c  is  a  constant 
=  134.68  ms;  e.g.  Wp  -  1503,  Bp  =4,  Ta  =  5865.72  ms  which  is  the  result  shown  in  Figure  2.6  for  Roll:  1, 
Xponden  1  Start:  14,  Proc#  2  in  the  Figure.  The  constant  c  is  a  sum  of  the  delays  through  the  array  system 
in  ms  determined  by  laboratory  tests  to  be  7  +  (128  -  0.016  *  (32  -  NPROC))  where  NPROC  is  the  number 
of  processors  in  the  array.  The  7  ms  is  due  to  data  buffering  in  the  array  hardware  and  the  remainder  is  due 
to  buffering  within  the  driver  for  the  magnetic  tape.  This  buffering  delay  was  discovered  later  to  be  in 
error  by  1 12  ps  and  should  have  been  calculated  as  (128  ms/frame  -  [125  (words/frame)  -  2  (header  words) 
-  NPROC  *  8  words/processor]  *  0.016  ms/word). 

The  result  of  the  correlation  is  seen  as  a  series  of  returns  across  the  array  for  each  of  the  transponders 
as  shown  in  Figure  2.7.  One  second  of  data  for  each  processes-  is  plotted  during  each  transponder  reply  with 
the  deepest  processor  PI  plotted  on  the  bottom.  The  reply  from  the  bottom  transponders  (Figure  2.7a,  b 
and  c)  appear  at  the  deepest  processor  first  and  arrive  sequentially  at  the  shallower  receivers  as  the  pulse 
travels  up  through  the  water  column.  The  pulse  from  the  last  transponder  (Figure  2.7d)  travels  down  from 
FLIP.  The  squared  amplitudes  are  normalized  to  1.0  and  plotted  such  that  a  value  of  1  will  be  slightly 
above  the  zero  level  of  the  next  processor.  The  noise  and  interference  mentioned  earlier  is  also  evident; 
processor  6  represents  a  particularly  noisy  time  series  and  imerferring  signals  are  seen  travelling  up  and 
down  the  array,  e.g.  between  0.6  and  0.7  s  on  the  plot  of  transponder  3.  Each  valid  data  window  is  marked 
with  ’Vs";  the  window  length  for  the  first  and  last  transponders  is  100  ms,  but  was  increased  for  tran¬ 
sponders  2  and  3  to  250  ms  to  enable  detection  of  the  transponder  error  modes  illustrated  by  the  green  tran¬ 
sponder  (transponder  2)  in  Figure  1.1.  The  time  of  arrival  of  the  detected  reply  may  be  estimated  from  the 
plot  by  adding  the  time  indicated  by  the  plot  to  the  number  of  ms  notated  at  the  top  of  the  plot  and  subtract¬ 
ing  a  constant  c  =  134  ms  defined  previously.  The  travel  times  are  converted  into  slant  ranges  by  assuming 
a  constant  sound  speed  of  1500  m/s  and  written  into  ascii  files  for  use  in  the  spatial  position  calculation 
computed  in  the  next  program. 


m.  Array  Spatial  Localization. 

The  navigation  algorithm  for  the  array  elements  is  virtually  identical  to  that  described  for  the  tran¬ 
sponders  once  the  initial  array  X-Y  position  is  determined,  except  that  the  transponders  are  considered  sta¬ 
tionary.  The  element  to  transponder  slant  ranges  corrected  for  sound  speed,  the  depth  of  the  transponder 
installed  on  FLIP  and  the  adjusted  transponder  locations  constitute  the  data  required  to  navigate  the  array. 
The  array  positions  are  iterated  to  achieve  the  best  fit  to  the  data  is  a  least  squares  sense.  A  description  of 
the  main  program  and  subroutines  and  a  sample  run  will  be  found  in  Appendix  C.  Array  element  relative 
location  accuracies  of  a  few  meters  are  achieved  by  this  method. 

The  array  is  not  moored  but  hangs  vertically  from  FLIP  under  320  kg  of  tension;  its  horizontal  range 
of  motion  is  significant  during  a  one  minute  time  interval.  Consequently,  receiver  locations  are  iterated 
using  a  single  time  slice  of  slant  range  data.  The  receiver  slant  range  is  the  path  from  FLIP  to  the  tran¬ 
sponder  to  the  receiver,  FLIP  slant  range  is  simply  the  path  to  the  transponder.  The  slant  ranges  are 
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TRANSPONDER  1 :  STARTING  AT  5893  ms  AFTER  ROLLOVER  BUFFER 


TRANSPONDER  2:  STARTING  AT  6048  ms  AFTER  BUFFER 


TRANSPONDER  3:  STARTING  AT  6138  ms  AFTER  BUFFER 


TRANSPONDER  4:  STARTING  AT  246  ms  AFTER  BUFFER 
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WINDOW 


Figure  2.7  Transponder  Reply  Detection.  The  correlator  output  for  each  processor  is  plotted  for  1  second 
around  the  time  of  each  of  the  transponder  replies.  The  top  processor  plotted  is  the  top  processor  in  the  array.  The 
valid  data  window  is  indicated  by  the  small  x’s.  The  bottom  mounted  transponder  replies  are  evident  as  upward 
traveling  pulses,  arriving  at  the  bottom  of  the  array  first  (a)  T1  (red)  transponder  reply  arriving  between  0.15  and 
0.55  for  processors  1  and  1 1  respectively;  (b)  T2  (green)  transponder  reply;  (c)  T3  (blue)  transponder  reply;  (d)  the 
12  kHz  pulse  transmitted  from  FLIP  arrives  as  a  downward  traveling  pulse  with  the  surface  bounce  120  ms  behind. 


21 


corrected  for  deviations  due  to  a  sound  speed  profile  which  differs  from  the  assumed  1500  m/s,  as  are  the 
receiver  depths  (which  are  obtained  from  the  travel  time  measurement  of  the  FLIP  transceiver  to  each 
receiver),  using  the  harmonic  mean  described  earlier.  The  noise  in  the  receiver  slant  range  measurements 
(Figure  3.1)  is  contributed  to  not  only  by  the  receiver  noise  as  shown  in  the  previous  section  but  by  the 
transponder  malfunction  described  in  the  first  section  and  by  high  frequency  FLIP  motion  (recall  that  FLIP 
slant  ranges  are  available  only  once  per  hour).  Constraints  on  the  difference  in  receiver  depths  and  slant 
range  measurements  from  one  minute  to  the  next  are  implemented  as  user  parameters  thresd  and  thres.  If 
these  thresholds  are  exceeded,  the  noisy  data  is  ignored  and  another  user  parameter  alter  determines  which 
interpolation  scheme  will  be  implemented,  if  any. 

Localization  of  the  array  receivers  now  proceeds  to  the  least  squares  filter.  A  constant  initial  xy 
position  is  assigned  to  FLIP  based  on  a  GPS  position  as  was  done  during  the  transponder  iteration  for  the 
FLIP  initial  position;  the  initial  positions  for  the  first  receiver  iterated  is  the  estimated  FLIP  position;  the 
initial  positions  for  subsequent  receivers  is  the  estimated  position  of  the  previous  receiver.  The  FLIP  slant 
range  is  subtracted  from  the  receiver  slant  ranges  leaving  the  transponder  to  receiver  portion.  The  horizon¬ 
tal  slant  ranges  are  calculated  and  the  xy  positions  are  iterated  minimizing  the  squared  error  in  calculated 
and  ’measured’  positions.  With  the  transponder  positions  fixed,  the  iteration  is  confined  to  FLIP  and  the 
array  receiver  positions  whose  x,  y  and  z  positions  are  output  with  an  associated  rms  error.  The  resulting 
analysis  of  data  processed  in  this  manner  is  presented  in  [Sotirin  and  Hildebrand,  1989]. 
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Figure  3.1  Slant  Range  Data.  The  slant  range  data  for  each  of  the  bottom  transponders  (as  detected  by  the  naviga¬ 
tion  receiver  sampled  by  processor  8)  are  plotted  relative  to  the  median  value  for  a  series  of  nonconsecutive  tapes. 
The  data  were  clipped  at  ±200  m  from  the  median  value.  The  errors  shown  in  Figure  1.1a  for  the  green  transponder 
are  evident  as  a  shadow  return  160  m  from  the  direct  return.  Most  of  the  data  appearing  at  -200  m  illustrates  a  lack 
of  detection  due  to  the  navigation  of  FLIP  once  an  hour  for  several  minutes.  FLJP  navigation  continues  to  interro¬ 
gate  the  transponders  consequently  the  12  kHz  returns  are  there,  but  because  the  FLIP  navigation  system  is  not  syn¬ 
chronous  with  the  array,  the  valid  data  windows  filter  out  most  of  the  unwanted  returns. 
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APPENDIX  A.  Transponder  Localization  Program  Description. 


The  tables  below  illustrate  the  program  flow;  the  text  details  purpose  of  each  routine.  The  structure  of  the 
tables  reflects  the  main  program  and  subroutine  levels.  Each  subroutine  in  level  1  is  called  by  the  main 
program;  each  subroutine  in  level  2  is  called  by  the  level  1  subroutine  to  the  left  The  structure  of  the  text 
is  similar,  indentations  reflect  subroutine  level.  Any  file  names  beginning  with  xxx  may  be  specified  by  the 
user.  The  first  routine  localizes  the  transponders  and  is  called  XPMA1N.  It  is  assisted  by  an  initilization 
program  called  XPLOAD.  These  programs  were  originally  written  for  navigating  the  transponders  used  in 
the  DEEPTOW  program  at  the  Marine  Physical  Laboratory.  A  sample  run  for  each  program  specifing  the 
input  and  output  file  structures  and  program  user  inputs  follows  the  description  of  the  program. 


XPLOAD 


Main 

Subroutines 

Program 

Level  1 

Level  2 

XPLOAD 

XPCMLD 

BASLIN 

XPSET 

XPCMLD 

EXIT 

XPLOAD:  initializes  the  program  parameters  used  by  the  main  transponder  navigation  program  XPMAIN. 
There  are  three  files  involved: 

XPCOMDAT  defines  a  common  area  containing  initialized  variables.  If  this  file  does  not  exist  when 
the  program  is  initiated,  it  is  created. 

xxx.trs  contains  the  X,  Y  and  depth  information  for  each  transponder  and  provides  a  convenient  way 
of  inputting  the  data.  The  default  name  is  TRANSPONDER.TRS,  read  in  if  a  carriage  return  is  input 
when  the  user  is  quenied  for  the  transponder  file. 

xxx  .1st  contains  the  X,  Y  depth  and  baseline  information  for  each  transponder.  The  user  is  querried  for 
the  name,  referring  to  the  ’listing’  file. 

XPCMLD:  is  a  subroutine  which  reads  or  writes  to  a  file  depending  upon  the  state  of 
Idflg.  The  file  is  called  XPCOMDAT  by  default  and  if  it  does  not  exist  when 
XPCMLD  is  called,  it  is  created. 

XPSET:  If  xpcom.dat  does  not  exist,  the  file  is  created;  program 
functions  include  initializing  the  deep  (1500  m/s)  and  shallow 
(1500  m/s)  sound  velocities,  logical  unit  numbers,  plotter  inputs 
and  clearing  the  transponder  data  buffer  areas.  Some  of  this  infor¬ 
mation  is  obsolete  eg.  sound  velocities  and  plotter  inputs. 

If  the  transponder  data  file  xxx.trs  does  not  exits,  it  is  created  by  quering  the  us~'.  The  inputs  for  each  tran¬ 
sponder  are  ’Label,  X,  Y,  Depth,  Comment’.  The  X  and  Y  values  are  positions  v.th  respect  to  an  arbitrary 
origin.  All  numeric  inputs  are  in  meters. 

BASLIN:  The  baselines  between  transponders  are  calculated  using  the  X-Y  positions 
and  are  written  into  the  list  file  along  with  the  X-Y  positions. 
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XPMAIN 


Main 

Subroutines 

Program 

Level  1 

Level  2 

XPMAIN 

XPCMLD 

XPINPT 

XXCOR 

LOOP: 

XPREAD 

XXCOR 

END  LOOP: 

XPLOOP 

XPURGE 

BASLN2 

TMDATE 

BASLN2 

BASLN2 

XPRINT 

TMDATE 

XPRIN2 

EXIT 

_ 

XPFIL  (perturbs  FLIP  position) 

XPFEL  (perturbs  transponder  positions) 

XPMAIN  perturbs  the  xy  positions  of  FLIP  and  each  transponder  until  the  RMS  error  is  minimized.  There 
are  four  files  involved: 

XPCOMDAT:  defines  a  common  area  which  is  set  up  by  XPLOAD.  The  common  area  contains  the 
transponder  information,  nominal  sound  speed,  plotter  parameters,  and  logical  device  definitions. 
Sound  speed  and  plotter  parameters  are  obsolete. 

xxx.dat:  contains  the  time,  X  and  Y  position  and  the  depth  of  the  interrogation  hydrophone  which  is 
deployed  from  FUP  and  slant  range  information  for  each  transponder  for  a  series  of  fixec.  This  file  is 
set  up  manually  by  reading  the  slant  range  for  each  transponder  off  of  the  chart  recorder  output  (if  sur¬ 
face  ship  is  used  for  the  survey,  by  also  recording  the  corresponding  LORAN-C  or  SATNAV  position 
and  converting  to  meters  from  an  arbitrary  origin).  The  same  xy  FLIP  position  is  used  for  all  fixes 
due  to  the  limited  range  of  the  survey. 

xxx.trsout  is  an  output  file  containing  the  adjusted  transponder  positions. 
xxx.pos  is  an  output  file  containing  the  adjusted  FLIP  positions. 

XPCMLD:  is  a  subroutine  which  reads  the  data  in  the  transponder  common  area  from 
the  file  called  XPCOMDAT  which  is  written  by  XPLOAD. 

XPINPT:  opens  the  data  file  xxx.dat,  the  transponder  file  xxx.trsout  and  inquires  as  to 
the  minimum  number  of  ranges  acceptable  (default=3).  The  option  of  selecting 
specific  transponders  to  be  removed  from  the  net  is  available  but  will  not  be  used  in 
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this  demonstration  as  the  minimum  number  of  transponders  is  3. 

XPMAIN  requests  an  operator  input  for  the  RMS  error  factor  for  fix  rejection  (default=2)  and  saves  the  ini¬ 
tial  transponder  positions  for  later  comparison. 

LOOP: 

XPREAD:  opens  the  data  file  xxx.dat  and  calculates  the  horizontal  range  (XY  projec¬ 
tion)  from  the  source  hydrophone  on  FLIP  to  each  transponder  and  for  each  fix  using 
the  source  depth,  transponder  depth  and  slant  range. 

HH=CRANS  (ntr  jipos  )=Vs2  -  D2 

where  ntr  indicates  a  particular  transponder,  npos  indicates  a  particular  fix  or  position, 

S  is  the  slant  range  from  the  source  to  that  transponder  and  D  is  the  transponder  depth  - 
source  depth. 

XPLOOP:  initializes  the  minimization  of  the  RMS  error  in  the  horizontal  range.  This  is 
accomplished  in  a  loop  in  which  XPLOOP  calls  XPFTL  transfering  the  appropriate 
parameters  to  first  adjust  the  FLIP  positions  (npos  fixes)  for  each  transponder  and  then 
adjust  the  transponder  positions  for  each  FLIP  position  (fix).  The  RMS  error  is  the 
error  after  perturbing  the  transponder  positions.  If  the  reduction  in  normalized  RMS 
error  from  one  loop  to  the  next  is  less  than  0.015  or  is  the  number  of  loops  exceeds  30, 
XPLOOP  is  terminated. 

XPFTL:  adjusts  the  xy  positions.  The  first  sequence  during  which  it 
is  called  within  the  XPLOOP  loop,  it  sums  the  squared  error  in  the 
xy  projections,  '£(rngec  -  HH )2  where  mgec  is  the  xy  projection 
calculated  from  the  xy  positions  of  the  FLIP  and  the  transponder 
and  HH  is  calculated  as  described  in  XPREAD,  for  each  tran¬ 
sponder.  The  FLIP  position  is  then  adjusted: 

X  =X  +DX  *  RATIO  *  GAIN  IN  where  X  is  the  original  X 
position  of  the  FLIP,  DX  =  the  difference  in  X-positions  of  the 
FLIP  and  the  transponder,  RATIO  *  (HH-RNGEQ/RNGEC, 
GAIN  =  1.5,  N  =  number  of  transponders.  If  the  root-mean  of  the 
squared  error  is  less  than  0.25  or  if  the  normalized  reduction  in  rms 
error  is  less  than  0.035  it  returns  to  XPLOOP.  Otherwise  the 
squared  error  is  cleared  and  calculated  again.  This  continues  until 
the  rms  error  satisfies  one  of  the  above  criteria  or  the  number  of 
adjustments  exceeds  30.  Although  the  final  error  calculated  from 
the  FLIP  position  adjustment  .s  passed  back  to  XPLOOP  it  is  not 
used.  The  parameters  associated  with  the  next  FLIP  position  are 
amassed  and  XPFIL  is  called  again.  This  continues  for  all  FLIP 
positions  or  fixes  in  xxx.dat.  The  second  sequence  during  which 
XPFIL  is  called  the  proceedure  is  essentially  identical  except  the 
the  squared  error  is  summed  for  each  FLIP  position,  the  position 
adjusted  is  the  transponder  position,  and  N  is  the  number  of  FLIP 
positions. 

XPMAIN  then  evaluates  the  rms  error  passed  back  by  XPLOOP  for  the  following  criteria: 

If  the  RMS  error  *  ^number  of  fixes  in  xxx.dat  is  <  1.0. 

If  the  normalized  reduction  in  RMS  error  per  loop  [(OLD  error  -  current  error)  /  OLD  error]  <  0.035. 

This  terminates  the  loop  if  the  reduction  in  error  is  less  than  a  specified  value  or  if  the  error  increases. 

If  the  RMS  error  is  less  than  a  specified  number  (0.75). 

If  any  of  the  above  criteria  are  met,  XPMAIN  jumps  out  of  the  LOOP  and  begins  the  print  sequence.  The 
error  must  converge  to  one  of  these  criteria,  there  is  no  limit  as  to  the  number  of  times  the  LOOP  may  be 
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executed. 


XPURGE  removes  any  FLIP  fixes  which  have  an  rms  error  greater  than  the  fix  rejec¬ 
tion  parameter  input  from  the  keyboard  earlier  (default=2)  times  the  total  transponder 
rms  error. 

BASLN2  calculates  baselines  from  transponder  xy  positions  and  writes  to  the  screen. 
XPMAIN  stores  the  new  transponder  positions. 

TMDATE  calculate  real  time  and  date  from  a  system  call. 

XPMAIN  writes  the  date,  time  and  transponder  positions  to  an  intermediate  file  which  is  defined  currendy 
as  the  screen. 

GOTO  LOOP 


BASLN2  calculates  baselines  from  initial  transponder  xy  positions  and  writes  to  the 
screen. 

BASLN2  calculates  baselines  from  final  transponder  xy  positions  and  writes  to  the 
screen. 

XPMAIN  calculates  the  difference  between  the  initial  transponder  baselines  and  the  final  transponder  base¬ 
lines  and  prints  to  the  screen. 

XPRINT  writes  to  the  screen  the  adjusted  transponder  ranges  and  the  time,  X,  Y,  error 
and  depth  for  the  FLIP  transceiver. 

XPRIN2  writes  to  the  output  file  xxx.pos  the  time,  X,  Y,  error  and  depth  for  the  FUP 
source,  a  quality  factor,  and  the  original  slant  ranges  to  the  transponders  in  the  same 
format  as  the  input  file  xxx.dat. 
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APPENDIX  B.  Array  Travel  Time  Acquisition  Program  Description. 


HARRYNAV 


Main 

Subroutines 

Program 

Level  1 

Level  2 

Level  3 

HARRYNAV 

getnav 

whichtp 

printnavp 

skipread 

gettime 

skipread 

dtfpos 

skipread 

fillnav 

dtfpos 

skipread 

printnavt 

whichtp 

navloc 

HARRYNAV:  Enters  user  parameters  from  the  program  line.  Parameters  are  the  option  flags,  number  of 
rollovers  to  skip  between  processing,  total  number  of  rollovers  to  process,  rollover  to  start  processing  with, 
processor  ID  file,  FLIP  slant  range  data  file  name,  and  output  file  name.  System  subroutines  are  called  to 
read  the  input  program  line  and  the  processor  ID  file. 

getnav:  System  subroutines  are  called  to  read  the  FLIP  slant  range  data,  parameters 
are  initialized.  Subroutines  whichtp  and  skipread  are  called  to  read 
in  the  buffer  of  data  containing  the  rollover  requested.  If  the  time  flag  is  set, 
the  first  buffer  on  the  tape  is  read,  the  header  time  parameter  is  converted  to  an 
integer  and  compared  to  the  user  input  time  parameter  also  converted  to  an 
integer  (by  gettime).  Due  to  a  lack  of  space  in  the  tape  header,  only  the 
least  significant  digit  of  the  GMT  hour  was  recorded  so  the  local  hour  is  con¬ 
verted  to  GMT  for  the  correct  GMT  time.  The  time  in  each  buffer  header  is 
compared  until  the  header  time  in  the  rollover  buffer  is  greater  than  the  time 
requested.  If  the  time  flag  is  not  set,  buffers  are  read  in  until  the  rollover 
number  requested  is  equal  to  the  number  of  rollovers  read  in.  A  loop  indexed  by 
i  over  the  number  of  rollovers  requested  is  initialized.  The  file  window. dat 
containing  the  valid  window  parameters  is  read  in.  Dtfpos  is  called  to  write 
the  date  and  time  into  the  output  file  and  FLIP  slant  range  into  the  output 
file.  A  loop  for  extracting  the  data  for  each  transponder  is  initialized  indexed 
byj.  Skipread  is  called  for  the  remaining  three  transponders.  Navtstl 
is  called  to  locate  the  return.  Files  are  output  containing  the  valid  data  markers 
( pltwinsio )  and  the  received  data  ( navplt.sio )  for  plotting.  Current  valid  data 
markers  are  written  into  window.dat  and  the  date  and  time  are  appended  to  the 
bottom  of  the  file  by  dtfpos.  A  system  call  to  pltsav.scr  plots  the  data  to 
the  screen.  If  pltsav.scr  does  not  exist  when  called,  a  message  is  displayed  and 
the  program  continues. 

whichtp:  Requests  the  user  to  load  the  magnetic  tape  and  enter  the  tape 
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drive  number.  If  the  verbose  flag  is  set,  printnavp  and 
printnavt  are  called. 

printnavp:  prints  the  array  slant  ranges  by  processor,  printnavt: 
prints  the  array  slant  ranges  by  transponder. 

skipread:  Reads  the  first  buffer  requested.  Buffers  may  be  requested  by 
rollover,  by  GMT  time  or  by  hardware  clock  time.  If  end  of  tape 
indicators  are  true  then  whichtp  is  called  and  the  program 
continues  extracting  the  navigation  time  series. 

gettime:  Translates  the  ascii  string  indicating  GMT  time  into  an  integer. 

dtfpos:  is  called  to  write  the  date  and  time  into  the  output  file  and 
linearly  interpolate  the  FLIP  slant  range  data  to  obtain 
the  slant  range  corresponding  to  the  time  of  the  rollover.  The 
local  date  is  converted  to  Julian  day,  and  the  time  is  GMT.  If 
the  passed  parameter  ’pall’  is  negative  (0)  then  only  the  data  and 
time  are  printed. 

fillnav:  For  each  transponder  interrogation,  the  navigation  bits  are 
extracted  from  the  data  stream  and  stored  as  a  time  series  for 
each  processor.  Each  time  series  begins  with  the  buffer  following 
the  interrogation  transmission  and  ends  74  buffers  later  (9.472  s). 
Various  ID’s  and  counters  are  verified.  Errors  in  the  buffer 
counter,  frame  sync  word  or  frame  counter  are  fatal.  If  an  error  in 
processor  ID  is  detected,  the  data  for  the  processor  in  question  is 
zero  filled. 

navloc:  The  data  are  received  by  navloc  as  the  most  significant  10  bits 
in  a  16-bit  word  and  repacked  into  the  30  least  significant  bits 
in  a  32-bit  word.  The  correlation  is  done  as  a  moving  adder  with 
leading  pointers  il  (word)  and  jl  (bit)  and  trailing  pointers  k 
(word)  and  1  (bit)  starting  from  the  beginning  of  the  plot 
window  (bwinw[xpndr]  +  bwinb[xpndr]  the  window  word  and  bit 
pointers)  for  SAVSZ  bits  (1  second  of  data).  The  plot  window  is 
set  up  as  WINOFF  words  before  the  valid  data  window  when  the 
program  is  initiated  and  allowed  to  move  as  the  reply  moves  by 
evaluating  the  average  movement  of  the  replies  across  the  array. 

The  detection  is  valid  only  within  the  valid  data  window,  with 
parameters  pstart,  pend,  LENGTH;  and  is  defined  as  a  normalized 
correlation  amplitude  of  detec=l  (normalized  by  PINGSIZE).  The 
time  of  the  detection  is  the  start  of  the  plot  window  plus  the 
number  of  bits,  p,  as  the  adder  was  initilized  with  the  first 
PINGSIZE  outside  the  correlation  loop. 
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406  4886.9  4890.3  4896.4  4900.7  4902.9 

411  4902.8  4903.0  4900.0  4902.1  4903.5 

416  4909.1  4909.1  4904.9  4901.4  4896.9 

421  4899.4  4897.1  4902.1  4898.6 

Sio  Channel  Number  4 


APPENDIX  C.  Array  and  FLIP  Localization  Program  Description. 


FLPNAV 

Most  of  the  code  for  this  program  is  involved  in  reading  the  inputs.  There  are  some  minor  system  subrou¬ 
tines  called  which  interpret  ascii  input  files.  Only  subroutines  which  pertain  to  the  array  navigation  prob¬ 
lem  will  be  described  here. 


Main 

Subroutines 

Program 

Level  1 

Level  2 

FLPNAV 

linint 

xcor 

fix3 

FLPNAV:  Locates  the  array  receiver  positions  in  a  least  squares  sense  using  parameters  entered  by  the  user 
in  a  command  file.  User  parameters  define  the  transponder  locations,  slant  range  input  data  file  name,  out¬ 
put  file  name,  sound  speed  profile,  interpolation  method  for  the  FLIP  slant  ranges  and  a  method  for  con¬ 
straining  deviations  in  slant  range  with  time.  The  set  of  measurements  taken  at  any  one  particular  time  are 
referred  to  as  a  "fix".  Some  parameters  may  be  input  from  the  command  file  or  from  separate  existing  files 
e.g.  dvfile  filename  inputs  depth  vs.  sound  speed  values  from  the  named  file.  Parameters  which  are  not 
specifically  declared  in  the  command  file  are  preset  to  the  value  indicated  below.  Some  of  the  parameters 
defined  below  were  used  for  test  purposes  only. 

The  information  required  by  this  program  is: 

1)  Fixed  x,  y,  and  z  positions  for  each  of  3  transponders. 

2)  Initial  x  y  positions  for  each  receiver  and  FLIP. 

3)  Fixed  depth  of  each  receiver  and  FLIP. 

4)  The  slant  range  of  each  receiver  and  FLIP  to  each  transponder. 

The  array  receiver  slant  ranges  must  be  the  distance 

from  FLIP  to  transponder  to  array  receiver. 

5)  Sound  speed  profile. 

The  navigation  procedure  is: 

1)  Discard  noisy  slant  ranges  using  the  parameters  THRESD  and  THRES. 

An  entire  fix  is  thrown  away  if  more  than  7  of  the  slant  ranges  for 

any  one  transponder  are  zero  (this  is  to  find  the  fixes  that  occur 
when  the  navigation  is  turned  off  -  usually  around  the  hour  when 
FLIP  is  navigated). 

2)  Adjust  the  receiver  depths  and  slant  range  measurements  using  the 
harmonic  mean  based  on  the  sound  speed  profile  measured  by  a  CTD. 

3)  Interpolate  for  a  new  FLIP  slant  range  set  if  files  are  available; 
subtract  the  FLIP  slant  range  from  the  receiver  slant  range  leaving 
the  range  from  transponder  to  receiver. 

4)  Calculate  a  horizontal  range  for  FLIP  and  each  anay  receiver 

5)  Iterate  the  FLIP  and  receiver  positions  minimizing  the  rms  error 
in  a  least  squares  sense. 
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Input  Parameter  Definitions: 

ifile  -  The  file  name  of  the  file  containing  the  slant  ranges  and  depths  of  FLIP 

and  the  array  receivers  whose  positions  are  to  be  iterated.  A  FLIP  slant  range 
is  the  path  from  FLIP  to  the  transponder,  a  receiver  slant  range  is  the  path 
from  FLIP  to  the  transponder  to  the  receiver. 

Preset  =  mglog.dat  e.g.  ifile  /localdata/ranges/data 
ofile  -  The  output  file  name.  When  given,  the  output  xyz  information  will  be  written  to 
file  ofile  instead  of  standard  out 
Preset  =  ’  ’  e.g.  ofile  output 

trlocs  •  A  list  of  the  transponder  locations,  given  as  x,  y,  z,  and  variance  for  each 
transponder. 

Preset  =  none  e.g.  trlocs  1000  1000  2000  1.0  3000  1000  2100  2.0  2000  2000  1500  1.0 
depths  -  A  list  of  depths  of  all  the  receivers,  replacing  the  receiver  depths  in  ifile. 

All  depths  including  FLIP'S  must  be  given  if  any  are  given.  DEPTHS  is  ignored 
unless  given. 

Preset  =  0  e.g.  90  300  375  450  525  600  675  750  825  900  975  .... 
nsta  -  The  number  of  fixes  in  ifile. 

Preset=  1 

stainc  -  The  increment  between  fixes  in  ifile  to  interate.  The  first  fix  in  ifile 
is  always  the  first.  e.g.  stainc  2  skips  every  other  fix. 

Preset  =1.  e.g  stainc  10 

n2ave  -  The  number  of  consecutive  fixes  to  average  before  trying  to  calculate  each 

location.  The  average  is  performed  by  doing  a  running  average  of  one-way  slant  ranges 
(transponder  to  receiver)  of  n2ave  fixes.  Noisy  fixes  (where  more  than  1/2  of 
the  ranges  are  "noisy")  are  no:  included  in  the  average.  n2ave  may  not  exceed  10. 
Preset  =  1  e.g.  n2ave  10 

alter  -  The  alternatives  for  noisy  slant  range  data.  Slant  ranges  more  than  thres 

different  from  the  slant  range  for  the  same  receiver/  transponder  of  the  previous 
day/time  will  be  "noisy". 

=0,  No  data  are  discarded  as  being  "noisy". 

=1,  A  single  noisy  slant  range  out  of  the  three  slant  ranges  for  a  single  receiver 
will  cause  the  data  for  that  receiver  to  be  discarded.  The  output  file  will  not 
contain  a  position  for  the  receiver  for  that  day/time.  The  file  will  still  contain 
the  total  number  of  receivers,  but  some  receiver  positions  will  be  duplicated  (this 
allows  plot  programs  to  just  plot  the  same  receiver  on  top  of  itself). 

=2,  The  noisy  slant  range  will  be  replaced  by  the  slant  range  for  the  same  receiver 
and  transponder  from  the  previous  day/time  in  the  input  file. 

=3,  The  "noisy"  slant  ranges  will  be  replaced  by  interpolating  the  slant  range  to  the 
same  transponder  from  adjacent  receivers.  A  noisy  receiver  on  either  end  of  the 
array  will  receive  the  slant  range  from  the  previous  fix.  Likewise,  if  two  adjacent 
receivers  have  noisy  slant  ranges,  the  first  one  will  receive  the  slant  range  for 
that  receiver  from  the  previous  day/time. 

=4,  Noisy  slant  ranges  are  calculated  by  extrapolation  or  interpolation  of  adjacent 
non-zero  slant  ranges. 

Presets  3 

thres  -  The  threshold  slant  range  in  meters  used  in  determining  noisy  slant  ranges. 

Preset  =  100.  e.g.  thres  50 

thresd  -  The  threshold  depth  difference  in  meters  used  in  determining  "noisy"  depths. 

A  depth  is  considered  noisy  if  it  is  more  than  THRESD  away  from  the  depth  of 
the  same  receiver  on  the  previous  day/time. 
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Preset  =  5 

trfile  -  A  file  containing  the  transponder  locations  ( x,  y,  z  )  in  the  format  written  by 
program  XPMAIN  with  an  additional  column  specifing  the  variance. 

Preset  =  none  e.g.  trfile  trs 

dvfile  -  A  file  containing  the  depth-sound  speed  pairs  as  described  below. 
vdp  -  A  list  of  sound  speed  depth-pairs  to  use  in  making  sound  speed  corrections.  The 
sound  speed  correction  is  made  by  first  convening  every  depth  and  slant  range  to  time 
using  a  constant  sound  speed  of  oldcv.  The  speeds  in  vdp  are  assumed  to  be 
the  speed  of  sound  at  the  depth  given.  The  speed  used  between  specified  depths  is  the 
"interval"  speed  which  is  the  average  speed  of  the  two  adjacent  depths. 
dvp  -  A  list  of  depth-speed  pairs.  This  is  identical  to  vdp,  except  that  the 
order  of  the  pairs  is  depth  followed  by  speed. 

Preset  =  none 

oldcv  -  The  constant  speed  used  in  obtaining  the  slant  ranges. 

Preset=  1500.  e.g.  oldcv  1450 

calctd  -  Assuming  that  the  array  is  vertical  and  directly  below  FLIP,  data  from  two 

array  receivers  may  be  used  to  estimate  the  transponder  depths.  Enter  the  two  array 
receiver  numbers  for  the  calculation.  Test  mode  only. 

Preset  =  0  e.g.  calctd  3  4 

Ipfile  -  The  pathname  of  an  output  file  containing  the  slant  ranges  in  the  format  expected 

by  the  looping  program  XPMAIN.  The  file  contains  the  slant  ranges  from  the  transponder 
to  each  receiver. 

Preset  =  none  e.g.  Ipfile  slr.loop 

erfile-  The  pathname  of  a  file  in  which  data  error  messages  will  be  written.  No  data  error 
message  will  be  written  unless  erfile  is  given. 

Preset  =  none  e.g.  erfile  oops 
fudge  -  A  factor  to  test  FLIP'S  slant  ranges.  Test  mode  only. 

Preset  =  0.  e.g.  fudge  10. 


Subroutines: 

linint:  Interopolates  the  FLIP  slant  ranges  linearly. 

xcor:  The  harmonic  mean  of  the  sound  speed  profile  is  used  in  correcting  the  slant 

ranges  and  depths  which  were  calculated  assuming  a  constant  sound  speed  of  1500  m/s. 
trgss:  Not  used  for  the  September  87  data.  Calls  Fix2  to  calculate  two 
possible  ’fix’  positions.  The  position  closest  to  the  third  transponder 
position  is  selected. 

fix2:  Nor  used  for  the  September  87  data.  Calculates  the  two  positions 
defined  by  the  intersection  of  the  slant  ranges  from  the  first 
two  transponders.  See  Appendix  C. 

fix3:  Adjusts  the  of  first  the  FLIP  position  and  then  each  array 

position  using  a  single  fix  to  obtain  a  least  squares  convergence. 

Outputs  the  adjusted  xy  position  and  the  rms  error  for  each  call. 
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119 . 10  1492.7  138.94  1487.5  158.78  1486.0  178.62  1485.3 
198.46  1484.9  218.30  1484.1  238.13  1483.7  277.79  1482.2 
297.62  1481.7  330.00  1481.2  350.00  1481.5  376.91  1481.3 
396.71  1480.8  416.54  1481.1  426.45  1481.2  456.17  1480.6 
495.79  1480.2  505.69  1479.4  525.50  1479.3  555.21  1479.2 


APPENDIX  D.  Initial  Position  Calculation. 


To  calculate  an  array  X-Y  position  (Figure  D  1 ),  the  horizontal  range  from  FLIP  or  the  array  to  any  tran¬ 
sponder  is  given  by 

Hproj  1  =  5  2  -  ( Dti-Dfup )2 

Hproj  1  is  the  horizontal  range  from  transponder  1  to  FLIP, 

S  is  the  one  way  slant  range  from  FLIP  to  transponder  1, 

Dt i  is  the  depth  of  transponder  1  and 

Dfup  is  the  depth  of  the  transmitter  mounted  on  the  bottom  of  FLIP. 

Two  X-Y  positions  are  calculated  from  the  intersection  of  the  arcs  defined  by  the  horizontal  projections: 
Xn,  Yj\  =  the  X-Y  position  of  transponder  1; 

Xt2,Yt2  =  the  X-Y  position  of  transponder  2; 

X  i2  and  T 12  are  the  differences  between  the  X-Y  transponder  positions  (X i2  =  XTX-XTi)\ 
Bline  =  the  transponder  baseline  range  =  (Xft  +  Y  fe)'A; 

Hproj  1  =  transponder  1  to  fix  range; 

Hproj2  =  transponder  2  to  fix  range; 

PI  =  projection  of  fix  onto  baseline  point  to  transponder  1 

P2  =  projection  of  fix  onto  baseline  point  to  transponder  2 

C  =  perpendicular  distance  from  fix  to  baseline 

XPj  -  X -coordinate  of  the  projection  of  the  fix  onto  the  baseline 

Ypj  =  Y-coordinate  of  the  projection  of  the  fix  onto  the  baseline 

Xl.Yl  =  X-Y  position  of  the  first  fix 

X2.Y2  =  X-Y  position  of  the  second  fix 


Bline  =Pl+P2 

C2  =  Hproj  \2  -  P  \2  =  Hproj!2  -P 22  =  Hproj!2  -  (Bline  -  P  l)2 


Hproj  l2-Pl2  =  Hproj!2  -  Bline2  +  2 P 1  Bline  -PI2 

p  j  _  Hpro  j  1 2  -  Hpro i !2  + Bline2 
IB  line 


COS0  = 


X 12  . 


sin8  = 


Y 12 
Bline 


Xpj  =Xt  i  +  P  lcos  0;  Ypj  —  Yt\  +  P  Isin  0 
X\=XPj-  Csin 0;  Y 1  =  Ypj  +  CcosQ 
X!  =  Xpj  +  Csin 0;  Y 2  =  Ypj  —  Ceos 0 
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Y 


Figure  D.l  Array  Position  Calculation.  The  parameters  required  to  calculate  the  lateral  position 
of  an  array  element.  The  points  at  FIX1  and  FIX2  give  the  intersection  of  the  horizontal  projec¬ 
tion  arcs  for  two  transponders.  A  third  transponder  range  is  used  to  choose  between  these  two  po¬ 
sitions. 
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APPENDIX  E.  Software  Listings. 


The  software  discussed  in  this  report  is  listed  below.  The  listings  contain  the  main  program, 
subroutines,  makefile  and  include  files  for  each  program  discussed. 
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eep  transponder  navigation  looping  program 


52 


53 


DO  110  NTR  -  l/PTTRS  CO  1830  I-2,NTRS 

•Save  "true"  positions  WRITE(LERR,822)  LABEDf I ) ,TRX(I ) ,TRY< I ) , 

IF  ( DUMMY ( 2 )  .eq.  1.0)  then  *  1*8(1) , (BSI£I£(N+J),J-1, 

call  xxcor(TRPN(NTR,3).0.0,1500.  ,0. 0, 0,TRPN(NTO, 3) ,  N  -  N  +  I  -1 

TRPN(WTR,3))  1830  CONTINUE 
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55 


56 


57 


59 


PRIWT  *,  ITIME,TIM(NPOS),<PSNS(NPOS,I),I-l,2),  *  ')  LESS  THAN  TRANSPONDER  DEPTH  (\F8.2, 

ERROR, DEPS, Q,  (RANS(WTR)  ,WTR-1,NTRS)  CDI11  FORMAT(16<  IX,  F14 . 2)  ) 


60 


61 


SUBROUTINE  XPRIN2 

XPRINT  PRINTS  OUT  TIME,  X,X,  AND  ERROR  TOR  FIXES 
INVOLVED  IN  TOE  LOOPING  OPERATION. 


XPINPT  (NAMFIE,IF1LE,EEILE,ITIMES) 


66 


IIMOVE  IS  ONE  WHEN  POSITION  IS  TO  BE  CHANGED. 


tjndate(datstr,  timstr) 


67 


ATDF  fortran  callable  version  of  'C'  routine  atof 
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(NdXNIW'(e)HVdSl)'(3H0W0N'  (T)HVdSI) 
(aaasKd'(T)3HVds) 


BJECTS  -  harrynav.o  navloc.o  fillnav.o  navtst2.o  reads. o  prints. 


lARRYNAV(l)  UNKNOWN  SECTION  OF  TOE  MANUAL  HARRYNAV(l) 
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Due  to  the  desire  to  keep  the  program  simple,  all  possible 
erroneous  input  conditions  have  not  been  guarded  against. 
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73 


74 


*  t d-trtiichtp( vf lg, nout , nroll , NPRQC , procnos , navch , *fd , t p) ; 
if  ( read (*fd, (char  *)u.cbuf ,bufsiz)  <  bufsiz) 

{  prlntf ("Fatal  tape  read  error  cntbuf-  %d\n",cntbuf ),-  exit(l); 
3wab( (char  *)u.cbuf,  (char  *)u.cbuf ,4 ) ,•  cntbuf++; 
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iprintf (fp,“  %10.3f“  ,c«a(i+3«Wr)+cl*a(i-l+3*Mr)); 

{print  f  ( fp,  "  %10. 3f  "  ,c«a  li*4*OTl +cl*a  (i-U4*VT) )  ,- 
c-89.92,  fprintf(fp,"  %10 . 3f\n“ ,c) ;  fflush(fp);  return 
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ERR  -  ERR/COUNT 

l average  error  squared  of  fix 

if  error  is  small  enough  or  is  no  longer  converging. ..quit 
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SUBROOTire  LINIOT(XA(YA,N,X,Y) 
DIMENSION  XA(N) , YA(N) 
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DIMENSION  vdp(ovdps) 
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